Quasiprobability distribution functions for finite-dimensional discrete phase spaces: 

spin tunneling effects in a toy model 
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We show how quasiprobability distribution functions defined over ./V 2 -dimensional discrete phase 
spaces can be used to treat physical systems described by a finite space of states which exhibit spin 
tunneling effects. This particular approach is then applied to the Lipkin-Meshkov-Glick model in 
order to obtain the time evolution of the discrete Husimi function, and as a by-product the energy 
gap for a symmetric combination of ground and first excited states. Moreover, we also show how an 
angle-based potential approach can be efficiently employed to explain qualitatively certain features 
of the energy gap in terms of a spin tunneling. Entropy functionals are also discussed in this context. 
Such results reinforce not only the formalism per se but also the possibility of some future potential 
applications in other branches of physics. 
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I. INTRODUCTION 



In the last decades, much effort has been devoted to 
characterize quantum tunneling processes in mesosco- 
pic and/or macroscopic systems, emphasizing the impor- 
tance of the degree of freedom related to the angular 
momentum and angle pair [l| . In what concerns the spin 
tunneling, certain theoretical approaches have pointed 
to some different ways of treating this problem, each one 
presenting a particular set of convenient inherent mathe- 
matical properties 0. From this perspective, if one con- 
siders physical systems with a finite-dimensional space 
of states and described by discrete variables, a sound 
theoretical framework must be employed to characterize 
properly such nonclassical effect. In fact, an alterna- 
tive approach to the system description of these spe- 
cific cases can be pointed out. First, we recognize that 
the state spaces associated with those particular physi- 
cal systems are TV-dimensional Hilbert spaces. Next, in 
connection with these finite Hilbert spaces, it should be 
stressed that quantum representations of _/V 2 -dimensional 
discrete phase spaces can also be constructed [3|. Thus, 
relevant operators whose kinematical and/or dynamical 
contents carry all the necessary information for describ- 
ing those quantum systems can now be promptly mapped 
in such phase spaces. In this sense, although there are 
various phase spaces formalisms proposed in the litera- 
ture for treating finite-dimensional physical systems 
let us focus our attention upon the framework devel- 
oped in Refs. 0, H, 0, H i] for the discrete represen- 
tatives of the quasiprobability distribution functions de- 
fined in 7V 2 -dimensional phase spaces, which has its al- 
gebraic structure based on the technique of constructing 



'Electronic address: mamarchi@ift.unesp.br 
^Electronic address: mandrake@ift.unesp.br 
t Electronic address: galetti@ift.unesp.br 



unitary operator bases initially formulated by Schwinger 
[icl | . The virtue of this discrete quantum phase-space ap- 
proach is that it allows us to exhibit and handle the pair 
of complementary variables related to a particular degree 
of freedom we are dealing with, as well as to recognize 
the quantum correlations between them. The basic idea 
then consists in exploring this mathematical tool in or- 
der to study those quantum correlations in connection 
with spin tunneling processes. Therefore, besides having 
the energy spectrum, which can be obtained via direct 
diagonalization of the Hamiltonian system, we can also 
have additional quantum information about the physical 
system through the study of the corresponding discrete 
Wigncr and/or Husimi functions. 

Such theoretical framework is then applied, in par- 
ticular, to the Lipkin-Meshkov-Glick (LMG) model pl| . 
which was originally introduced over forty years ago in 
nuclear physics [12| for treating certain fermionic sys- 
tems. This important toy model has, since then, been 
extensively studied in the literature because of its appa- 
rent simplicity [l3[ . Indeed, it can also be viewed as a fi- 
nite set of spins half mutually interacting in the xy-plane 
subjected to a transverse magnetic field [l4j]- Moreover, 
our interest in this model also resides in the fact that 
spin tunneling can be considered to occur. 

In this work, we show how the time-dependent dis- 
crete Husimi function and an angle-based potential des- 
cription [TEl . [l6l . IT?! ] can be combined in order to describe 
a group of physical processes that encompasses, among 
other things, the spin tunneling effects. This particular 
approach leads us not only to extract the energy gap for a 
symmetric combination of ground and first excited states, 
but also to corroborate its inherent applicability to analo- 
gous physical systems such as in magnetic molecules [18| . 

This paper is organized as follows. In Section II, we 
present a condensed review of the theoretical apparatus 
used to describe 7V 2 -dimensional phase spaces and also 
discuss some new essential features exhibited by the time- 
dependent discrete Husimi function. In Section III, we 
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apply our results to the LMG model to explore qualita- 
tively the spin tunneling effects for a symmetric combi- 
nation of ground and first excited states. Finally, Section 
IV contains our summary and conclusions. 



II. THEORETICAL APPARATUS FOR 
FINITE-DIMENSIONAL PHASE SPACES 

Our theoretical framework is totally based upon the 
formalism developed in Refs. 0, [|[ for physical systems 
with finite-dimensional space of states. In this sense, let 
us introduce some basic elements which represent our 
guidelines for the fundamentals of the formal descrip- 
tion of iV 2 -dimensional phase spaces by means of discrete 
variables. The first important element is the mod(TV)- 
invariant operator basis 



riA=- 



which consists of a discrete Fourier transform of the ex- 
tended mapping kernel S( s '(r),£) = [J^(rj, £,)]~ s S(r], £), 



Sfa,0 := -L expfee j U"V : 



(2) 



being the symmetrized version of the unitary operator 
basis proposed by Schwinger [1(3]. In this particular ap- 
proach, the labels 77 and £ are associated with the dual 
momentum- and coordinate-like variables of a discrete 
7V 2 -dimensional phase space. Note that these labels are 
congruent modulo N and assume integer values in the 
symmetrical interval [—(,£) for £ = (N — l)/2 fixed. Be- 
sides, the extra term J^(r], £) is defined through the ratio 
M{r\, £)/^#(0, 0), where the function Q 

+ ^ 4 (ary|ia) [i9 3 (ae|ia)e i7T « + tf 4 (a<e|ia)e i7r (" + « +Ar >] } 

with a = (27V) -1 , is responsible for the sum of products 
of Jacobi theta functions evaluated at integer arguments, 
and s refers to a complex parameter that satisfies the re- 
lation \s\ < 1. It is worth mentioning that a compilation 
of results and properties which characterize the algebraic 
structure of the discrete mapping kernel T( s ) (fi, v), as 
well as the unitary operators U and V, can be promptly 
found in Refs. p, 0|- For physical applications related 
to quantum tomography and quantum teleportation, see 
alsoRef. 0. 

Next, let us assume that pit) reflects the dynamics of 
a particular quantum system characterized by a finite- 
dimensional space of states. The one-to-one mapping 
between density operators and functions belonging to 
an A^-dimensional phase space labelled by {/z, v} £ 
[—£, £] is attained, in our description, by means of the 



parametrized function (fj,,v;t) := Tr[T' s )(^, v)p{t)}. 
Expressed as a double discrete Fourier transform of the 
discrete s-ordered characteristic function S( s ) (77, £;t) := 
Tr[S( s )(?7,£)p(r;)], it has a well-established continuous 
counterpart within the Cahill-Glauber formalism 
Moreover, for s = —1,0,1 the time-dependent func- 
tion [jj,, v\ t) is directly related to the respective dis- 
crete Husimi, Wigner, and Glauber-Sudarshan distribu- 
tion functions. An interesting formal result from this for- 
malism explores the connection between discrete Wigner 
and Husimi functions - here denoted by W(/x, v\t) and 
!K(/U, v ; t) - through the equation 



M(A*,i/;t) = l E( M ,H//y)W(//y ; i), (3) 

where ~&{p,,v\p! ,1/) = Tr [Tl^/j^jT^^'y)] defines 
a smoothing process characterized by a discrete phase- 
space function that closely resembles the role of a Gaus- 
sian function in the continuous phase space. In fact, 
Eq. ([3]) represents an intermediate smoothing sequence 
within a hierarchical process among the quasiprobability 
distribution functions in finite-dimensional spaces |8[ . In 
what concerns the time evolution of the density operator, 
it must be stressed that p{t) satisfies the von Neumann- 
Liouville equation and its corresponding mapped expres- 
sion in the discrete phase-space representation leads us to 
obtain a differential equation for W(u,^;i), whose solu- 
tion was explicitly determined in Ref. [f|. Hence, the dis- 
crete Husimi function can now be immediately inferred. 

Now, let us establish an important mathematical re- 
sult for the time-dependent discrete Husimi functions. 
First of all, it should be noticed that JC(/x, v\ t) is strictly 
positive and limited to the interval [0, 1] for any t > 0; 
secondly, the phase space treated here consists of a finite 
mesh with iV 2 points and characterized by the discrete 
variables [i and v. So, the discrete Husimi function can 
be mapped onto a real NxN matrix H(t) whose elements 
[h rs (t)] r s=1 N obey two essential properties that imply 
in the conservation of probabilities on a discrete phase 
space, that is, 

N 

< h rs (t) < 1 and ^ h rs(t) = 1 • 

r,s=l 

The interaction with any dissipative environment is au- 
tomatically discarded within this context. Such mathe- 
matical procedure brings some operational advantages in 
our description since the matrix H(t) can be promptly 
diagonalized onto the eigenspaces 

V Al = {¥,; eV:HV i = AiVi} 

characterized by the eigenvectors {V^t)}^!,...^ an d 
their corresponding eigenvalues {\i(t)}i=i t ...,N for t fixed. 
It is worth emphasizing that the eigenvalues obtained 
from this particular diagonalization process assume, in 
general, both real and complex values, and this fact can 
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be explained by means of matrix analysis [2(j. A perti- 
nent question then emerges from our considerations on 
7V 2 -dimensional phase spaces: "Can both real and com- 
plex eigenvalues be associated with some physical pro- 
cess?" 

To answer this question, let us initially decompose the 
matrix H(t) as a sum of two Hermitian and antihcr- 
mitian N x N matrices, namely, M(t) = A(t) + B(i), 
both matrices being constructed out following the ma- 
thematical recipe A(i) = \ [H(f) +H f (t)] and B(i) = 
| [H(i) — Ht(t)l. In this situation, the diagonalization 
process attributes real eigenvalues for A(t), while B(i) 
has eigenvalues which are pure imaginary (or zero). Be- 
sides, the trace of H(i) is preserved, i.e., 

N 

Tr[H(t)]=Tr[A(t)]=^<7«(t). 

»=i 

Here, {a^ (t)}i=i,...,jv and {fg (t)}»=i,...,jv represent the 
respective eigenvalues of the matrices A(t) and B(i) for 
all t > 0; therefore, it is easy to verify that these eigenva- 
lues are now responsible for the real and imaginary parts 
of the complex eigenvalues {Ai(i) }*=!,. ..,jv. Next, let us 
introduce an auxiliary tool characterized by the entropy 
functional 

N 

E[{A i };i] = -^|A i (t)|ln(|A i (t)|) , (4) 

i=l 

which allows us to infer all the different contributions as- 
sociated with {Xi(t)} £ C. From the operational point 
of view, the simplicity of this measure represents an ef- 
fective gain to our task since the diagonalization process 
of the time-dependent discrete Husimi function is suffi- 
cient in this case for determining E[{Aj};f]. The next 
step then consists in considering a well-known physical 
system that leads us to find out any concrete evidence 
in ([3]) of some particular physical process inherent to the 
model. In this sense, we will apply the theoretical frame- 
work here discussed to a solvable quasi-spin model whose 
Hamiltonian, although simple, presents some interesting 
physical and mathematical features, namely, the Lipkin- 
Meshkov-Glick (LMG) model [TO]. 

III. THE LMG MODEL 

Originally proposed with the intent of testing mean- 
field approximations in many-body systems, the LMG 
model is here introduced through the Hamiltonian [n| 

H = 2 X] ^ a l°- a 1> CT + ~2 a io- a q',o- a l''-°" ac l'- cr 

which describes a collection of N p fermions distributed 
in two Np-fold degenerate levels separated by an energy 
e. The degenerate states within each level are labelled in 



this expression by means of the quantum numbers q € 
[1,N P ] and a — ±1 (+1 and —1 represent the respective 
higher and lower levels), N p being considered an even 
number. 

It is worth noticing that the introduction of the quasi- 
spin operators |12| 

J± := a l±i a q<Ti and J 3 : = \ CT a 3,,<T a q,^ 

q q,<T 

into the LMG Hamiltonian not only reveals its underly- 
ing SU(2) structure, since the operators J± and J 3 obey 
the standard commutation relations [J+,J_] = 2 J , and 
[J 3 , J±] = ±J±, but also allows us to treat collective ex- 
citations of the fermionic system in a more suitable form, 
namely H = eJ } + (V/2) (j\ +J 2 _). Indeed, the 3 } term 
of this Hamiltonian operator gives half the difference of 
the number of particles laid on the upper and lower le- 
vels, while the second term, involving the operators 
and J_, is associated with the interaction between a pair 
of particles located at the same energy level, being also 
responsible for the scattering process of this pair to the 
other level where the quantum number q of each particle 
is preserved. For simplicity, let us rewrite such Hamil- 
tonian operator in order to scale the interaction term to 
the particle number while measuring the energy in terms 
of e, i.e., 

H ' : =7 = J ' + 2l;« +j2 -> < 5 > 

with x = N p y/e. Since [Hi, J 2 ] = 0, it turns immediate 
to see that Hi can be diagonalized within each (2J + 1)- 
dimensional multiplet labelled by the eigenvalues of J 2 
and J 3 , which accounts for the soluble character of the 
associated quantum model [lj, [T3| ■ It is worth stressing 
that the ground state belongs to the finite multiplet cha- 
racterized by J = = max(J 3 ), so that its J 2 and J 3 

quantum numbers are (-^j flf + -*-) and — if > respec- 
tively. Hereafter, we will be interested only in this par- 
ticular multiplet containing the ground state and whose 
dimension is given by N = 2J + 1 = N p +1. This (N p + 1)- 
dimensional multiplet will be considered as our underly- 
ing Hilbert space of interest [21| , and also will be used to 
construct the (N p + l) 2 -dimensional discrete phase space. 

Next, let us mention some few words about two dis- 
crete conserved quantities inherent to the LMG model 
which reflect certain symmetry properties. The simplest 
operator commuting with the Hamiltonian (O, there- 
fore giving a constant of motion, is the parity operator 
n := exp(i7rJ 3 ). This fact tells us that the Hamiltonian 
matrix, in the J 3 representation, breaks into two disjoint 
blocks involving only even and odd eigenvalues of J 3 , res- 
pectively. The second interesting quantity comes from 
the anticommutation relation {Hi,R} = 0, where 

R(— 7r/2, 7r, 0) := exp[i(7r/2)J 3 ] exp(i7rJ !/ ) 

corresponds to a rotation of the angular momentum 
quantization frame by the Euler angles (— tt/2, it, 0), thus 
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(a) 




20 



FIG. 1: Time evolution of "K(m, n; r) for the Lipkin-Meshkov- 
Glick model with N p = 20 and \ = 1.5 fixed, where the labels 
m and n characterize, respectively, the dimensionless angular 
momentum and angle pair. These pictures show, in particu- 
lar, how the two-body interaction term present in the Hamil- 
tonian H affects the initial distribution "K(m,n;0) for diffe- 
rent values of dimensionless time t. We have adopted in our 
numerical investigations the values (a) r = 0, (b) r « 6.5, (c) 
r « 15.9, and (d) r « 25.3, which illustrate a representative 
but not complete evolution of the discrete Husimi function. 



transforming — > Hl- In this case, if |Ej) is an 
energy eigenstate with eigenvalue E J ; then R|Ej) is also 
an eigenstate of with eigenvalue — Ey. This symme- 
try property of the Hamiltonian operator © gives rise 
to an energy spectrum that is symmetric about zero. 

After this condensed review, let us establish a se- 
quence of steps that permits us to evaluate the time 



evolution of the discrete Husimi function for the LMG 
model. The first one consists in adopting the theoretical 
approach developed in Ref. [l5| for the time-dependent 
discrete Wigner function W(m, n; t) defined upon an N 2 - 
dimensional phase space labelled by the angular momen- 
tum and angle pair (m,n) G [— N p /2,N p /2]. Since this 
approach depends on the Wigner function evaluated at 
time t = 0, the second one consists in fixing the ini- 
tial state as a symmetric combination of the ground and 
first excited states following the recipe described in 
The next and last step refers to the smoothing process 
given by Eq. ^ , which leads us to finally obtain the de- 
sired result. Figure 1 shows the three-dimensional plots 
of J{(m,rt;r) versus (m, n) € [—10,10] with N p = 20, and 
X = 1.5 fixed. In the numerical investigations, we have 
adopted some specific values for the dimensionless time r 
in order to illustrate the effects of the two-body interac- 
tion term in the original Hamiltonian H (or, cquivalently, 
the second term of constituted by the operators H\ 
and J"l) on a given initial configuration of the finite phase 
space. Thus, figure 1(a) represents the discrete Husimi 
function !H(m, n; 0) which reflects the initial condition of 
the model under investigation, that is, it shows two dis- 
tinct regions equally distributed in the 21 2 -dimensional 
phase space of the particular two-level system. Figure 
1(b) corresponds to a subsequent time r « 6.5, where 
we perceive that the probability distribution was almost 
totally reallocated in one side of the discrete phase space, 
which means that both the angular momentum and angle 
components present negative values. By its turn, figure 
1(c) shows an intermediate process for r w 15.9 and quite 
similar to that found in (a) when r = 0. Finally, figure 
1(d) illustrates the migration, in r w 25.3, to a distribu- 
tion of positive values of angle components in constrast to 
that viewed in figure 1(b). Note that, in particular, the 
quantum dynamics of this system inhibits a distribution 
function centered at n — 0. 

Although the periodic pattern verified for the discrete 
Husimi function can be explained, in principle, through 
the periodic fluctuation of particle populations between 
the two energy levels, it is not clear until now its relation 
with the energy gap A = (Ei — Eo)/e, associated with the 
quasi-spin tunneling occuring in the present situation, as 
well as its dependence on the interaction parameter \- To 
clarify this point, some considerations on the energy spec- 
trum related to deserve be properly mentioned: (i) it 
is numerically obtained by just diagonalizing the matrix 
associated with the Hamiltonian (O in the 3 l -basis (for 
more details on the exact solutions of the LMG model, 
see Ref. [HI); as a direct consequence of this result, 
(ii) the energy gap curve shows an explicit dependence 
on x, that is, as the two-body correlation strength in- 
creases, the energy gap decreases [Til, 13 ■ Therefore, if 
one changes the interaction parameter \ ( or the energy 
gap A), the periodic pattern showed by J{(m,n;r) will 
be also modified. Next, with the help of the entropy func- 
tional ([3]) , we will estimate this quantity for the same set 
of parameters used in the previous figure. In this way, 
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FIG. 2: Time evolution of E[{A;};t] versus r G [0,60] for 
the same set of parameters used in the previous figure. It is 
worth noticing that the maximum and minimum points of this 
curve correspond to the cases when the discrete Husimi func- 
tion reaches its peaks in the respective negative and positive 
regions of the angle sector of the finite phase space. This fact 
permits us to estimate the energy gap A = (Ei — Eo)/e from 
our numerical data, namely A w 0.1784, which corroborates 
the value obtained through the diagonalization process of the 
Hamiltonian operator ((5} in the J,-basis (A = 0.1788). 
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Figure 2 shows the plot of E [{Ai};r] versus r £ [0,60] 
for N p = 20 and \ — 1.5 fixed. Note that the entropy 
functional curve presents an oscillatory behaviour with 
a well-defined periodic structure, which allows us to es- 
timate the energy gap through the period of oscillation 
between two consecutive maximum (minimum) points. 
Thus, after a detailed analysis of the numerical data used 
in the plot of Fig. 2, we can conclude that A « 0.1784. 
It is important emphasizing that this result is in perfect 
agreement with that obtained from the diagonalization 
process for (in this case, the percent error estimated 
is <5 w 0.22%). 

Finally, we will present some plausible arguments that 
lead us to establish, under certain circumstances, a link 
between the oscillations of the discrete Husimi function, 
the energy gap, and the spin tunneling effect. For this 
intent, let us initially adopt the theoretical framework 
exposed in Ref. [3J, where the potential function 



V{4>) 



1 



COS(0) + -: 



2N, 



■ sin 2 (0) 



(6) 



and the 'effective mass' function 



M-\4>) = ^— r + X [1 + sin 2 (0)] } (7) 

were derived with details and exhaustively tested for the 
LMG model in the special limit N p 3> 1. Such functions, 
defined in the interval cj> € [— tt, tt], permit us to explain, 
from a phenomenological point of view, the underlying 
behaviour of the discrete Husimi function observed in 
Figure 1. For instance, figure 3(a) shows the plots of 
V(<p) (solid line) and M{4>) (dashed line) versus <fr for 
N p = 20 and x = 1-5 fixed; in addition, figure 3(b) repre- 
sents an amplified image of the curve related to V (</>), 
where now the lowest energy levels E /e = —10.9343 
and Ei/e = —10.7555, extracted from the diagonaliza- 
tion process of the Hamiltonian in the J 3 -basis, are 



FIG. 3: Figure 3(a) represents the plots of V(cj>) (solid line) 
and M{4>) (dashed line) versus cf> £ [— with N p = 20 
and x ~ 1-5 fixed. Note that both curves have similar be- 
haviours and share complementary informations about im- 
portant quantum effects such as that appeared in figure 3(b), 
where now the energy levels Eo/e and Ei/e have also been 
drawn. The possibility of tunneling through the potential 
barrier can be blocked if the effective mass function presents 
a divergent behaviour at or near the classical turning point; 
otherwise, the tunneling effect takes place (see both pictures 
at 4> = 0). Such description of spin tunneling, via potential 
and effective mass functions, can also be extended to magnetic 
molecules, such as Mnl2-acetate molecule and Fe8 clusters. 



also viewed (see dot-dashed lines). It should be stressed 
that distinct values of x modify not only the energy gap 
but also the behaviour of the curves related to the po- 
tential and effective mass functions. Hence, the value 
here chosen for x brings out some explicit advantages 
of this theoretical approach since the potential function 
presents a pronounced barrier at the origin which affects, 
consequently, the energy levels Eo and Ei . Indeed, figure 
3(b) consists of a paradigmatic case where tunneling ef- 
fects take place in this quasi-spin system. Within this 
context, figures l(a,c) would then correspond to a spin 
wavepacket equally distributed in both sides of the sym- 
metric double- well potential, while figures l(b,d) repre- 
sent a wavepacket localized only on one side of the po- 
tential barrier centered at <f> = 0. Therefore, the time 
evolution of J£(m, n; r) can now also be understood as 
a faithful representation upon a discrete iV 2 -dimensional 
phase space (here labelled by the dimensionless angular 
momentum and angle pair) of spin tunneling processes 
associated with the symmetric combination of the ground 
and first excited states. 
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FIG. 4: Plot of l[Ji;r] versus r £ [0,60] for N p = 20 and 
X = 1.5 fixed. The oscillatory pattern viewed in this picture 
also allows us to estimate the energy gap from our numerical 
data, i.e., A » 0.1786, with an error of 5 » 0.11%. 

IV. CONCLUDING REMARKS 

In what concerns the symmetry property related to the 
energy spectrum of the LMG Hamiltonian, let us men- 
tion some few words about the lowest energy doublet. 
Numerical investigations have shown that, if one consi- 
ders the energy symmetric counterpart of this doublet, 
the time evolution of the spin wavepacket will reveal a 
pronounced oscillatory behaviour, and this result is di- 
rectly associated with the high energetic demand of the 
system in accessing these particular states. The smooth- 
ing process due to Eq. (J3]) and the subsequent diagonali- 
zation of the discrete Husimi function will then produce, 
for all r > 0, an entropy functional E[{Ai};r] with a 
well-defined periodic pattern but not related to the pre- 
viously discussed tunneling effects. Therefore, the theo- 
retical apparatus here exposed shows that the coherent 
oscillations verified in this case do not characterize the 
same spin tunneling process, although the value of the 
energy gap be the same. 

Next, we introduce a complementary functional to that 
established by E [{A^}; t], which permits us to measure, in 
principle, the correlation between the discrete variables 
of an iV 2 -dirnensional phase space. To this end, one con- 
siders the mutual correlation functional [t| H2] 

I[IK; t] := E[Q; t] + E[ft; t] - E[IK; t] > , (8) 

where E[IK;i] corresponds to the time-dependent joint 
entropy defined in terms of the discrete Husimi function 
(J3J), with E[Q; t] and E[IR; t) representing the marginal en- 
tropies which are related to the respective marginal dis- 
tributions Q(n;t) and "R{v;t) (for technical details, see 
Refs. [H, |||). Thus, if one applies this measure to the 
LMG model, some interesting results can be promptly 
obtained. In this sense, figure 4 illustrates the time evo- 
lution of I[5f;r] versus r e [0,60] for the same set of 
parameters fixed in the previous figures. From the nu- 
merical point of view, it is immediate to perceive that: 
(i) the maximum points coincide with the configurations 
exhibited in figures l(a,c) and also reflect a situation 
where the spin wavepacket is equally distributed in both 
sides of the potential barrier, which implies in a mini- 



mal mutual correlation (maximal uncertainty) between 
the angular momentum and angle pair; (ii) the minimum 
points describe both the configurations illustrated in fi- 
gures l(b,d) and this result can be explained by means of 
a spin wavepacket localized in one of the potential wells, 
which leads us to obtain a maximal mutual correlation 
(minimal uncertainty) between the discrete variables m 
and n; and finally, (iii) I [IK; r] presents a half period if 
one compares with the oscillatory pattern exhibited by 
E[{Ai}; t], this fact being associated with the absence of 
differentiation between the cases reported in l(b,d). 

In summary, we have developed an alternative theo- 
retical framework for a class of physical systems de- 
scribed by discrete variables with potential applications 
in quantum information theory and quantum computa- 
tion [23J . Based on a finite-dimensional phase space des- 
cription, this formalism was then applied to the Lipkin- 
Meshkov-Glick model [ll[ whose Hamiltonian operator, 
although apparently simple, presents some remarkable 
physical and mathematical properties 13]. In particular, 
we have shown how the angle-based potential approach 
US EE E3 can be used to explain qualitatively the spin 
tunneling effects related to a symmetric combination of 
ground and first excited states. Moreover, we have also 
inferred as a by-product the energy gap for this situation 
(i.e., related to the particular spin tunneling situation 
described here) through two different ways, and showed 
that both methods produce excellent quantitative results 
if one compares them with the exact analogues extracted 
from the diagonalization process of the LMG Hamilto- 
nian. 

Finally, it is worth mentioning that an interesting 
study about Fe8 magnetic clusters in the presence of ex- 
ternal magnetic fields appeared recently in [24| , where the 
theoretical apparatus for finite-dimensional phase spaces 
here discussed was extensively applied with great success 
in exploring the spin tunneling effects in those clusters (in 
particular, such quantum effects occur for temperatures 
below the crossover temperature, i.e., 0.35-ftf [25j|). Since 
the Fe8 magnetic clusters are characterized by a J = 10 
spin ground state [III, some similarities with the LMG 
model can be shared via discrete phase-space approach. 
However, it must be stressed that such phenomenological 
Hamiltonian, describing the Fe8 magnetic clusters, con- 
tains essential experimental information concerning the 
anisotropies inherent to the molecule structure [27j , what 
produces an energy spectrum that is essentially diferent 
from that obtained for the LMG model. From this theo- 
retical approach, it was also possible to infer the energy 
gap between the ground and first excited states, namely 
A sa 2.922411^, with a percent error 6 estimated in 1.04% 
(in this case, the external magnetic field applied along 
the easy axis of the cluster has as intensity Hn = 0.1 IT). 
Besides, the entropy functional E[IK; r] and I [IK; r] were 
also employed to qualitatively explain how the spin tun- 
neling effect is connected with the functional correlations 
between the discrete variables observed in the underlying 
phase space, corroborating, in this way, the results here 
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